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ABSTRACT 

We present a semi-analytic model for the shaping of the Homunculus Nebula around Carinae that 
accounts for the simultaneous production of bipolar lobes and an equatorial disk through a rotating 
surface explosion. Material is launched normal to the surface of an oblate rotating star with an initial 
kick velocity that scales approximately with the local escape speed. Thereafter, ejecta follow ballistic 
orbital trajectories, feeling only a central force corresponding to a radiatively reduced gravity. Our 
model is conceptually similar to the wind-compressed disk model of Bjorkman & Cassinelli, but we 
modify it to an explosion instead of a steady line-driven wind, we include a rotationally-distorted star, 
and we treat the dynamics somewhat differently. A continuum-driven explosion, where the radiation 
force is independent of velocity, avoids the disk inhibition mechanisms that normally operate in line- 
driven winds. This allows mid-latitude material with appropriate initial specific energy to migrate 
toward the equator where it collides with material from the opposite hemisphere to form a disk. 
Thus, our model provides a simple method by which rotating hot stars can simultaneously produce 
intrinsically bipolar and equatorial mass ejections, without relying on an aspherical environment or 
magnetic fields. Although motivated by rj Carinae, the model may have generic application to episodic 
mass ejection where rotation is important, including other luminous blue variables, B[e] stars, the 
nebula around SN1987A, or possibly even bipolar supernova explosions themselves. In cases where 
near-Eddington radiative driving is less influential, our model generalizes to produce bipolar pinched- 
waist morphologies without disks, as seen in many planetary nebulae. If rotating single stars can 
produce strongly axisymmetric ejecta by this mechanism, then the presence of aspherical ejecta by 
itself is insufficient justification to invoke close binary evolution. 

Subject headings: circumstellar matter — ISM: individual (Homunculus Nebula) — stars: individual 
{rj Carinae) — stars: mass loss — stars: rotation — stars: winds, outflows 



1. INTRODUCTION 

It is commonly assumed that shell nebulae surround- 
ing massive hot stars like luminous blue variables (LBVs) 
consist of slow ambient material that has been swept-up 
by the faster wind of the hot supergiant. This scenario is 
often adopted to explain the origin of bipolar geometry 
in their nebulae, by applying the generalized interacting 
stellar winds (GISW) scenario that was developed suc- 
cessfully for bipolar planetary nebulae and related out- 
flow phenomena (e.g., Mellema et al. 1991; Frank, Balick, 
& Livio 1996; Balick & Frank 2002), as weU as for the 
famous nebula around SN1987A (Luo & McCray 1991; 
Wang & Mazzati 1992; Blondin & Lundqvist 1993; Mar- 
tin & Arnett 1995). In this scenario, a hot fast wind 
expands into a slower wind from a previous red super- 
giant (RSG) or asymptotic giant branch (AGB) phase. 
The surrounding slow wind must have an equatorial den- 
sity enhancement (i.e. a disk), and the consequent mass 
loading near the equator slows the expansion of the shock 

^ Based in part on observations obtained at the Gemini Observa- 
tory, which is operated by AURA, under a cooperative agreement 
with the NSF on behalf of the Gemini partnership: the National 
Science Foundation (US), the Particle Physics and Astronomy Re- 
search Council (UK), the National Research Council (Canada), 
CONIC YT (Chile), the AustraUan Research Council (Australia), 
CNPq (Brazil), and CONICET (Argentina). 



interface between the two winds, giving rise to a pinched 
waist and bipolar structure. However, it remains unclear 
how the required pre-existing disk can be formed. One 
does not normally expect RSG or AGB stars to be ro- 
tating rapidly enough to form an equatorial decretion 
disk such as is characteristic to Be stars (Porter & Riv- 
inius 2003), and thus a disk-shedding scenario probably 
requires the tidal influence of a companion during prior 
evolutionary phases in order to add sufficient angular mo- 
mentum. In the case of SN 1987A, a binary merger would 
be required for this particular scenario to work (Collins 
et al. 1999; Morris & Podsiadlowski 2006). As an alter- 
native, variations of the GISW paradigm can be invoked 
(e.g., including an aspherical fast wind expanding into a 
spherical slow wind) to reproduce the bipolar shape of 
the nebula around rj Carinae (Frank et al. 1995, 1998; 
Dwarkadas & Balick 1998; Langer et al. 1999; Gonzalez 
et al. 2004b, 2004b). 

However, it is unlikely that this general scenario can 
work for massive LBVs like 77 Carinae, or for SN1987A. 
Stars with luminosities above roughly 10^'^ L0 never 
reach the RSG stage, and the coolest apparent tempera- 
tures that they can achieve occur instead during the LBV 
phase. ^ As they evolve off the main sequence, they move 

^ Admittedly, there is also a group of relatively low luminosity 
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to the right on the HR diagram, toward eooler temper- 
atures, larger stellar radii, and lower values for their es- 
cape velocities. Consequently, their stellar wind speeds 
get slower — not faster — as their mass-loss rates in- 
crease. 0-type stars have typical wind speeds of a few 
10^ km s~^, whereas LBVs typically have terminal wind 
speeds of a few 10^ km s^^ and mass- loss rates a fac- 
tor of ^100-1000 higher. This creates a situation where 
a slow dense wind is expanding freely into a faster and 
much lower-density wind, which is exactly the opposite 
situation of that required for the usual GISW scenario 
to work. In other words, the winds are not strongly in- 
teracting. 

The specific case of SN1987A presents its own set 
of difficulties, even though its progenitor is well below 
log(L/LQ)=5.8 and it probably has been through a re- 
cent RSG phase. First, a merger model followed by a 
transition from a RSG to BSG requires that these two 
events be synchronized with the supernova event itself 
(to within the ~10'' yr dynamical age of the nebula), re- 
quiring that the best observed supernova in history also 
happens to be a very rare event. One could easily argue, 
though, that the merger (needed for the bipolar geome- 
try) and the blue loop scenario might not have been in- 
vented if SN1987A had occurred in a much more distant 
galaxy where it would not have been so well-observed (i.e. 
we wouldn't know about the bipolar nebula or the BSG 
progenitor). Second, after the RSG swallowed a compan- 
ion star and then contracted to become a BSG, it should 
have been rotating at or near its critical velocity (e.g., 
Eriguchi et al. 1992). Even though pre-explosion spectra 
(Walborn et al. 1989) do not have sufficient resolution to 
measure line profiles, Sk-69°202 showed no evidence of 
rapid rotation (e.g., like a B[e] star spectrum). Third, 
and particularly troublesome, is that this merger and 
RSG/BSG transition would need to occur twice. Prom 
an analysis of light echoes for up to 16 yr after the su- 
pernova, Sugerman et al. (2005) have identified a much 
larger bipolar nebula with the same axis orientation as 
the more famous inner triple ring nebula. If a merger 
and RSG/BSG transition are to iDlame for the bipolar- 
ity in the triple-ring nebula, then what caused it in the 
older one? Perhaps a more natural explanation would 
be that Sk 69°202 suffered a few episodic mass ejections 
analogous to LBV eruptions in its BSG phase (see Smith 
2007). The B[e] star R4 in the Small Magellanic Cloud 
may offer a precedent at the same luminosity as the pro- 
genitor of SN 1987A; R4 is consistent with a 20 M© evolu- 
tionary track, and it experienced a minor LBV outburst 
in the late 1980's (Zickgraf et al. 1996). R4 also has 
elevated nitrogen abundances comparable to the nebula 
around SN 1987A. 

For rj Carinae, observations have falsified the idea that 
the bipolar shape arises from a prior equatorial density 
enhancement (Smith 2006). In the GISW scenario, the 
nebula's waist gets pinched because of mass loading, so 
in the resulting bipolar nebula there should be an excess 
of mass at low latitudes compared to a spherical shell. 
Instead, detailed observations of the latitudinal mass de- 
pendence in the Homunculus show that the mass was 

LBVs around log(L/Lo)=5.5 (see Smith, Vink, & de Koter 2004a) 
where the GISW may still apply, because these stars may be in a 
post-RSG phase. 



concentrated toward high polar latitudes (Smith 2006). 
Furthermore, observations show no evidence for a pre- 
existing slow disk around tj Car that could have pinched 
the waist. Instead, the disk structure that is seen in high 
resolution images, usually called the "equatorial skirt", 
has measured kinematics indicating an origin at the same 
time as the bipolar lobes (Morse et al. 2001). Some ma- 
terial appears to be even younger, not older (Smith & 
Gehrz 1998; Davidson et al. 1997, 2001; Smith et al. 
2004b; Borland et al. 2004). Thus, any model for the 
production of the bipolar Homunculus also needs to be 
able, simultaneously, to produce a thin equatorial disk. 
The only model proposed so far to accomplish this is 
in the thermally-driven magnetic wind model of Matt & 
Balick (2004) intended for the presently-observed bipolar 
wind (Smith et al. 2003a). However, it is not known if 
the extreme magnetic field required to shape the massive 
Homunculus is achievable, because the conditions during 
the 19th century eruption were much more extreme than 
in the stellar wind seen now. 

If LBV nebulae cannot be dominated by swept-up am- 
bient material, we are left in need of an alternative expla- 
nation for the origin of their bipolar structure. Several 
recent clues, due in large part to detailed observations 
of T] Carinae, point instead toward the idea that LBV 
eruptions can behave more like explosions than steady 
winds: 

1. The ratio of total mechanical energy to radiated en- 
ergy for the 19th century eruption of r] Car was greater 
than unity, and the ratio of kinetic momentum to photon 
momentum was ~10^ (Smith et al. 2003b). These num- 
bers are characteristic of explosions rather than winds. 

2. Most of the mass in the Homunculus is concentrated 
in a very thin shell seen in H2 emission (Smith 2006), 
which implies that At for the mass-loss event was <5 yr. 
Proper motions also suggest a similarly small range of 
ejection dates (Morse et al. 2001). Ejecting the observed 
amount of mass in the Homunculus over that short a time 
requires a mass-loss rate of order 1 M© yr"-*^ or higher 
(Smith et al. 2003b; Smith 2006), which may surpass the 
capability of even a super-Eddington continuum-driven 
wind by itself (Owocki et al. 2004). 

3. The high polar expansion speed in the Homunculus 
(about 650 km s^^) is close to the expected escape ve- 
locity of the primary star ignoring radiation pressure. In 
a radiatively-driven stellar wind, one expects the speed 
to drop as a star nears the Eddington limit because of 
the lower effective gravity. However, in an optically-thick 
hydrodynamic explosion that is not initially drivc^n by ra- 
diation force, this condition would not necessarily apply. 

In this paper we examine the idea that a surface explo- 
sion from a rapidly-rotating hot star can create an intrin- 
sically bipolar nebula and thin disk, without relying on 
any latitude dependence in the ambient material, binary 
influence, or magnetic fields. Thus, our model adopts 
the simplest set of assumptions that are also physically 
plausible. 

2. THE MODEL: A ROTATING SURFACE EXPLOSION 

2.1. Basic Principles 

Our simple model traces its origins back to a key ques- 
tion concerning the morphology of 77 Car's Homunculus: 
can a single basic physical paradigm account for the si- 
multaneous production of both the equatorial skirt and 
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Fig. 1. — A simple conceptual sketch of the various levels of complexity for trajectories in a rotating explosion, (a) Purely radial 
trajectories from a spherical non-rotating star, (b) Orbital trajectories from a spherical rotating star, where material is diverted toward 
the equator as in the wind compressed disk model of Bjorkman & Cassinelli (1993). (c) Same as Panel a, but with ejection normal to the 
surface of an oblate star, (d) A combination of Panels b and c, where material is ejected normal to the surface of an oblate star, but where 
rotation modifies the trajectories. This still enhances the polar density, but rapid rotation also diverts material toward the equator to form 
a compressed disk. The shading in Panels c and d is to remind us of equatorial gravity darkening. 



the bipolar lobes? In considering this question the phe- 
nomenon of rotation naturally springs to mind, since 
material kicked ballistically from the surface of a ro- 
tating star has a natural tendency to migrate toward 
the stellar equatorial plane. This notion was invoked in 
the wind-compressed disk (WCD) model of Bjorkman & 
Cassinelli (1993) as a means for explaining the circum- 
stellar disks around Be stars. The WCD model treats the 
wind plasma as independent fluid parcels whose dynam- 
ics are confined to their respective orbital planes. De- 
tailed calculations by Owocki, Cranmer & Gayley (1998) 
demonstrated that subtle line-driving effects associated 
with velocity-dependent forces tend to inhibit the for- 
mation of disks in the WCD model. However, in our 
treatment below we focus on continuum driving, which 
is independent of velocity, and the WCD model furnishes 
a useful paradigm on which to build. 

We advance the hypothesis that, to first order, the 
shape of the Homunculus can be understood as the conse- 
quence of the nearly ballistic, anisotropic ejection of the 
surface layers of a star near its critical rotation limit. In 
this scenario, the bipolar shape and equatorial disk arise 
as a direct consequence of the ejection geometry, rather 
than through hydrodynamic interaction afterward. We 
assume that the subsequent evolution of these layers is 
governed solely by a central, gravity-like (1/r^) force, 
and therefore that the trajectory followed by each fluid 
element may be treated as a problem in orbital dynam- 
ics. The exception is for elements that pass through the 
equatorial plane; there, we assume that they collide with 



material from the opposite hemisphere, leading to the 
cancellation of their velocity component normal to the 
plane. We ignore the self gravity of the ejected mass. 

Sketches of the trajectories for outflowing material un- 
der various levels of complexity are illustrated in Figure 
1. These trajectories show in principle how density can 
be enhanced at the equator to form a disk, as well as 
toward the poles. When a spherical star with radial tra- 
jectories (Fig. la) is allowed to rotate rapidly, material 
is diverted toward the equator (Fig. 16) by conservation 
of angular momentum, as envisioned in the WCD model. 
However, a severely oblate star introduces other effects as 
well. First, as noted above, an oblate star has a smaller 
radius, larger effective gravity, and larger escape veloc- 
ity at the pole than at the equator. This - rather than 
the latitude dependence of density ~ is what gives rise to 
the bipolar shape of the nebula. The oblateness of the 
star will also enhance the polar density if material is ini- 
tially launched normal to the surface, simply because the 
flatter poles of the star aim trajectories poleward (Fig. 
Ic). We do not attempt to model this enhanced polar 
density quantitatively, however, because we expect that 
additional effects may be as important. For example, the 
polar mass flux can be increased by the stronger radiative 
flux at the poles of a star that suffers from rotationally- 
induced gravity darkening (e.g., Owocki & Gayley 1997; 
Owocki et al. 1998). In any case, the straight (i.e. non- 
rotating) trajectories from a rotationally-distorted star 
in Figure Ic seem unrealistic unless the ejecta speeds are 
much faster than the rotation speed. The more likely 
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rotationally-modified trajectories are shown conceptually 
in Figure Id. This last case creates a compressed disk, 
but also retains some degree of enhanced density at the 
poles, because the rotational wind compression effects 
are more important at larger cylindrical radii. 

Elaborating on this sketch, we consider a model star 
in an initial state of critical rotation. (We have no di- 
rect evidence that 77 Car was rotating critically prior to 
the Great Eruption; however, since the investigation of 
rotational effects is central to our study, it is appropri- 
ate to focus on the critical limit in which these effects 
are the most pronounced.) In the Roche approximation 
(e.g., Cranmer 1996, and references therein), the surface 
radius of the oblate star is given by 



f o sin(9/3) n () ^ /o 
R(f)) = I sine ^ TT/Z 

^ ' \B.{tT-0) e>TT/2 



(1) 



where 6 is the usual polar coordinate, and Rp is the polar 
radius of the star. The azimuthal velocity due to rotation 
is given by 

V^{9)^ncR{0)sm0, (2) 



where 



nr. = 



l8GM 
27i?3' 



(3) 



is the critical angular frequency of rotation in the same 
Roche approximation, with M the stellar mass and G 
the gravitational constant. 

At t — 0, we disrupt this initial state by (i) impart- 
ing a velocity kick V\^{d) to each surface fluid element 
in the direction of the local surface normal (this is the 
explosion), and (ii) introducing a spherically symmetric 
force that everywhere is directed radially outward, and 
whose magnitude is a times the gravitational force. (It 
is assumed that a < 1, so that the net force on elements 
remains directed inward toward the star). This disrup- 
tion represents our basic characterization of the Great 
Eruption observed in the 1840's, when the nebular ma- 
terial is thought to have been launched (Morse et al. 
2001; Smith & Gehrz 1998; Currie et al. 1996; Gehrz & 
Ney 1972; Ringuelet 1958; Gaviola 1950). For simphcity 
we assume instantaneous ejection rather than sustained 
outflow over 10-20 years, but this should have little effect 
on the overall results. A series of bursts over a decade, 
as opposed to one single burst, would produce a nebula 
with a similar shape but some finite thickness. 

Our motivation for including the kick (i) and the sub- 
sequent outward force (ii) comes from considering the 
effects of continuum radiation driving during the sud- 
den, factor 5 increase in luminosity associated with the 
Great Eruption (see, e.g., Davidson & Humphreys 1997). 
At the beginning of the eruption, rotation-induced grav- 
ity darkening will produce a strongly anisotropic radia- 
tion field (e.g., von Zeipel 1924; Cranmer & Owocki 1995; 
Owocki et al. 1996, 1998; Owocki & Gayley 1997; Langer 
1998; Glatzel 1998; Maeder 1999; Maeder & Desjacques 
2001). This means that the additional radiative flux es- 
caping the star will at first deposit momentum prefer- 
entially at the stellar poles; by assuming an appropriate 
form for Vk{0) (discussed in greater detail in the follow- 
ing sections), we use the velocity kick (i) to model this 
initial polar deposition. As the surface fluid elements 



subsequently move outward, however, the anisotropics 
in the radiative flux will tend to be smeared out, leading 
to a more spherically-symmetric outward radiative force 
that we incorporate via (ii). 

To determine the t > evolution of the surface fluid 
elements, we assume that each follows a trajectory de- 
scribed by the equation of motion 



dv 

dt 



CM 



GM 



a- 



(4) 



Here, r is the position vector of the element, v = dr/dt 
the corresponding velocity vector, and we adopt the con- 
vention that non-bold symbols denote the modulus of 
their bold vector equivalents, so that in this case r = |r|. 
The acceleration terms on the right-hand side of this 
equation arise, respectively, from the inward gravita- 
tional force and the outward spherically-symmetric force 
introduced at t = 0. Following the discussion given 
above, the initial velocity of each element is calculated 
as 

vo = Vk(e)e„-Hy0(0)e^ (5) 

where 

(6) 



e, - {R'/R) ee 
VI + [R'/RY 



is the unit surface normal vector, with R' = dR/d6, and 
{er,eg,e^} are the unit basis vectors in the spherical- 
polar {r, 9, (j)} directions. 

The equation of motion ^ is identical to that for a test 
particle moving in the gravitational fleld of a point mass 
M{1 — a) > 0. Therefore, the solutions are analytical, 
taking the form of conic sections (ellipses, hyperbolae, 
etc.) whose focus lies at the stellar origin. We leave a 
detailed discussion of these solutions to any of the many 
standard texts discussing this classical two-body prob- 
lem (e.g., Boccaletti & Pucacco 1996). However, it is 
appropriate to specify how the six orbital parameters, 
defining the trajectory followed by each surface fluid el- 
ement, are determined from the initial conditions. First, 
we calculate the specific (per-unit-mass) angular momen- 
tum vector 

j = ro X vo, (7) 

(where Tq is the element's position at t = 0, and Vq was 
defined in cqn. and the specific Laplace- Runge-Lenz 
vector 

A = vx j-GA/(l-a)e^. (8) 

We assume that the reference plane is the stellar equato- 
rial (Cartesian x — y) plane, with the x-axis defining the 
Vernal point. Then, the inclination i of the orbital plane 
is given by 

j ■ ■ (9) 



cos I 



J 



the longitude of the ascending node Q and argument of 
periastron lu by 

j ■ [A X (e^ X j)] ^ 
jA ■ (e^ X j) ' 
(10) 

and the eccentricity e and semi-major axis a by 



tan CI = 



Gz ■ [ex X (e^ X j)] 
• (gz X j) 



tancj = 



GM{1 - a) ' 



GM{1 ~ a) 



- 2GM{1 - a)/r 



(11) 
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(In these expressions, {ej:,ez} are the unit basis vec- 
tors in the Cartesianja;, z} directions.) Finally, the true 
anomaly vq of the surface element at the t = epoch is 
given by 



The orbital parameters {i, fl, uj, e, a, vq\ defined above 
allow calculation of the complete i > evolution of 
a given surface element. However, for those elements 
whose trajectories pass through the equatorial plane, the 
parameters must be modified to account for the antic- 
ipated collision with material from the opposite hemi- 
sphere. As discussed above, the polar {9) component of 
the velocity v is set to zero when the element reaches the 
equator. Then, the orbital parameters are recalculated 
using the element's instantaneous position and updated 
velocity, and its evolution is continued. 

This approach is rather different from the WCD model 
of Bjorkman & Cassinelli (1993), who assumed that ve- 
locity vectors become radial at the equator, with no 
change in speed. In fact, a more fundamental differ- 
ence between our approach and the WCD model lies in 
the treatment of the outward force introduced at t = 0. 
Bjorkman & Cassinelli (1993) incorporated a parameter- 
ization of line-driven wind theory in their model, result- 
ing in an outward radiative force that (a) exceeds grav- 
ity, and (b) does not have a simple 1/r^ scaling. The 
significance of (a) is that an initial kick (as assumed in 
our treatment) is not required for material at the stel- 
lar poles to escape from the star. However, (b) means 
that the equation of motion does not correspond to a 
two-body gravitational problem, and must be integrated 
numerically. 

2.2. Illustrative Simulations 

To furnish an initial demonstration of our approach, 
we consider the case where the velocity kick function is 
given by 

H(f?) -H,o|cos0|, (13) 

for some normalizing velocity Vk,o &t the poles, which in 
the case of r] Car is about 650 km s~^ (Smith 2006). The 
dependence on | cos 6*1 may seem ad hoc in our simulation, 
because it essentially prescribes the overall shape that is 
observed. However, this latitude dependence for the ini- 
tial kick has a firm physical justification. As we have 
discussed above, gravity darkening — whereby 
the local emergent flux scales with the local ef- 
fective gravity — tends to initially focus the ad- 
ditional radiative flux escaping from the star to- 
ward the stellar poles. Thus, the kick imparted 
by this flux is expected to be strongest over the 
poles, suggesting the above form for Vk{9). 

We should mention that a prescription similar to (|13p 
above has already been invoked to explain the bipolar 
shape of the Homunculus Nebula and the latitude de- 
pendence of rj Car's stellar wind for near-critical rotation 
(Owocki 2005; see also Owocki & Gayley 1997; Maeder 
& Desjacques 2001; Dwarkadas & Owocki 2002; Smith 
2002, 2006; Smith et al. 2003a). However, the focus in 
most of these studies is on the terminal velocity Vaa of 
a line-driven wind, whereas in the present study we are 
considering the initial kick velocity of an explosive ejec- 
tion. The reason why a cos9 variation is appropriate in 



the line driven case is that Voo typically scales with the 
local escape velocity Vcsc (see, e.g., Dwarkadas & Owocki 
2002). Then, with Vesc itself scaling with effective grav- 
ity in the same way as the radiative flux, a coincidence 
between Vk and the line-driven Voo naturally arises. 

For a selected region of a-Vk,o parameter space, we 
conducted simulations where we evolve a set of surface 
elements to the asymptotic limit t — > oo. For each simu- 
lation, the initial state at i = is comprised of 5,000 ele- 
ments distributed uniformly in 9 over the Roche surface 
described by cqn. ([1]). With initial velocities described 
by eqns. (|2l5ll3p . these elements are evolved according to 
the equation of motion ([5]), as described in the preceding 
section. In the limit of large t, the envelope Re{9) defined 
by the elements reaches a steady state that corresponds 
to the eventual shape of the ejected nebula. 

Figure m compares simulated envelopes Re{9) against 
the shape of the ry Car Homunculus as measured by Smith 
(2006). Clearly, with an appropriate choice of parameters 
— in this case, a = 0.6 and Vk,o — 1-5 ^JGM/Rp — 
we are able to capture the gross qualitative features of 
the bipolar nebula, while at the same time producing 
the desired equatorial skirt. Generally, a skirt occurs 
whenever the initial specific energy of surface elements, 

^0 = |vo|V2-GM(l-a)/ro, (14) 

exhibits a minimum at some point between equator and 
pole. (Such minima themselves arise because the kick 
kinetic energy decreases toward the equator, while the 
rotational kinetic energy and gravitational potential en- 
ergy both increase toward the equator.) The surface ele- 
ments situated at the energy minimum lead to the narrow 
waist of the envelope. The elements closer to the equa- 
tor, with higher Eq and hence faster terminal velocities 
Woo oc \/E^, then form the skirt, while the elements closer 
to the poles produce the bipolar lobes. In Fig. [3] we plot 
both i?o and \/E^ as a function of for the a = 0.6, 
Vk_o — 1-5 \/GM/Rp simulation. The specific energy 
minima that give rise to the skirt in this case can clearly 
be seen at 6* « 80° and 9 « 100°. Note that the tendency 
of elements to evolve toward the equatorial plane, due to 
conservation of angular momentum, produces a rather 
narrower skirt than might be assumed from a cursory 
look at the plot; compare, for instance, the \/Eq data 
shown in Fig.[3l with the center panel of Fig. [2] Figure 3 
reproduces the trend of higher kinetic energy in the polar 
ejecta observed in the Homunculus (Smith 2006). 

Toward smaller values of a and/or Vk,o the initial spe- 
cific energy of surface elements near the equator is neg- 
ative. These elements therefore remain gravitationally 
bound to the star, producing a collapsed- waist morphol- 
ogy with no skirt. The implication is that, when gen- 
eralized to cases where near-Eddington radiative driving 
is less influential [a <0.5 in Fig. 2), our model simpli- 
fies and easily accomodates the more common pinched 
waist morphologies of planetary nebulae without fast 
skirts (e.g., Balick & Frank 2002). This is encouraging, 
since these sources have central stars that are indeed far 
less luminous than rj Carinae. Conversely, toward larger 
a the elements have insufficient angular momentum to 
reach the equatorial plane, and the skirt is replaced by a 
local inversion of the bipolar shape. This range of condi- 
tions may explain why equatorial disks like that around 
rj Car are not always seen, while bipolar ejecta nebulae 
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Fig. 2. — The envelopes He(6) of the ejected surface elements in the asymptotic limit t ^ oo (solid), plotted together with the measured 
shape of the Homunculus (dotted) from Smith (2006). Each panel shows results from a simulation with the indicated gravity reduction 
parameter a and velocity kick normalization Vkfi, the latter being measured in units of ^GM/ Rp. 



are quite common around evolved rotating stars. 

2.3. A Tuned Simulation 

A notable discrepancy between observations and our 
'best' (central panel) simulation shown in FigureOis that 
the apparent flattening of the nebula lobes over the poles 
is not properly reproduced. This can be seen as faster 
expansion in the Homunculus at latitudes about 15-40° 
from the polar axis, as compared to our predicted shape. 
It is as if the Homunculus has received an extra kick at 
these latitudes. 

What Vk(6') must be adopted in order to reproduce 
correctly the lobe shapes? One possible answer to this 
question is presented in Figured Here, the kick velocity 
function has been tuned with the specific purpose of re- 
producing the lobe shapes. The tuning was accomplished 
through a simple process of trial and error, subject to the 
constraints that a — 0.6 and that over the stellar poles 
Fk = 1.5 ^GM/Rp (th ese values come from the 'best' 
simulation of the preceding section) . The close match be- 



tween theory and observations seen in the central panel 
of Figure |4] should not be taken as a measure of the fi- 
delity of our model — indeed, since we have adopted an 
ad hoc prescription for Vk(^), the close agreement is to 
be expected. 

However, what is of particular significance is the fact 
that the simulation simultaneously reproduced the lobes 
and a skirt that is similar to that which is actually seen 
around 77 Car (e.g., Duschl et al. 1995). Namely, this is 
a true flattened disk-like structure, with material of the 
same age spread out over a range of radii in the equa- 
torial plane - it is not a ring, even though we adopted 
an instantaneous explosion event (it is also kinematically 
different from a Keplerian disk). This is a true success 
of our model, since the signature of this skirt is wholly 
absent from the initial kick velocity function in the left 
panel of Figure 3. 

Furthermore, the small extra kick needed at latitudes 
about 15-40° from the polar axis may have a reasonable 
physical explanation. Although our model attributes the 
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Fig. 3. — The initial specific (per-unit-mass) energy Eq (sofid), 
and its square root V Eg (dotted) , plotted as a function of 9 for the 
a = 0.6, Vtfi = 1.5 simulation displayed in Fig. [2] 

overall shape of the Homunculus and its disk to initial 
conditions of the ejection, this is a case where hydrody- 
namic shaping of the ejecta long after ejection may play 
some role after all. We have argued that the GISW sce- 
nario cannot drive the overall shape of the Homunculus, 
largely because the post-eruption wind is not powerful 
enough (Smith 2006; Smith et al. 2003a, b). However, the 
post-eruption wind may be able to modify or perturb the 
existing shape. It is often seen in hydrodynamic simula- 
tions of interacting winds that bipolar nebulae develop 
"corners" at the outer parts of the polar lobes, leading to 
flattening over the poles that is reminiscent of the extra 
velocity we require here (e.g., Cunningham et al. 2005; 
Frank et al. 1998; Dwarkadas & Bahck 1998). This ef- 
fect arises in the following way: at low latitudes within, 
say, 45° of the equator, the fast wind that is inflating the 
bipolar nebula strikes the inner side walls of the dense po- 
lar lobes at an oblique angle, and is then deflected pole- 
ward. This material skims along the inner side walls of 
the polar lobes and piles up in the corners (see Cunning- 
ham et al. 2005), adding an extra kick over a small range 
of latitudes. Although the average momentum of r] Car's 
post eruption wind is tiny compared to the momentum 
of the nebula (Smith et al. 2003b), in this scenario, the 
momentum from a relatively large volume fraction of the 
post-eruption wind is focused on only a small portion of 
the nebula. This effect is purely hydrodynamic, so our 
ad hoc approach in Figure 3 may not be the best way to 
explore it. A better (but computationally more intense) 
way might be to take our predicted model shape in the 
central panel of Figure 2, and allow it to be "inflated" 
and shaped self-consistently by a post-eruption wind. 

The potential role of the post-eruption stellar wind 
(which is a line-driven wind) that follows the continuum- 
driven outburst also suggests an interesting "double 
whammy" for enhancing the bipolar shapes of LBV neb- 
ulae. Specifically, Owocki and collaborators (Owocki & 
Gayley 1997; Owocki et al. 1996, 1998) have demon- 
strated that radiative line driving not only inhibits the 
initial formation of an equatorial disk, but also enhances 
the wind mass-loss rate toward the stellar poles. In the 
case of T] Car, this polar wind has been observed (Smith 



et al. 2003a) . Such a polar wind may enhance the bipolar 
shape of the initial ejection. 

Of course, one can imagine other factors that may 
modify the shape of the polar lobes in the required way, 
such as the influence of a nearby companion star (e.g., 
Morris & Podsiadlowski 2006). This realm of shaping 
mechanisms is beyond the scope of our present study, 
but should be explored further. 

3. DISCUSSION 

Using a semi-analytic model, we have shown that a 
surface explosion from an oblate star near critical rota- 
tion can simultaneously produce an equatorial disk and a 
pair of polar lobes that closely approximate the observed 
shape of the Homunculus Nebula around t] Carinae. This 
model is arguably the simplest model that also includes 
realistic assumptions. It does not require any effects of 
hydrodynamically interacting winds or magnetic fields 
to produce the asymmetry. This shows that rotating hot 
stars can eject intrinsically bipolar nebulae simultane- 
ously with equatorial disks. 

Our techniques combine two aspects of theories devel- 
oped initially for non-spherical line-driven stellar winds, 
but we modify them to the scenario of a sudden explo- 
sion. Namely, our model borrows conceptually from the 
WCD model of Bjorkman & Cassinelli (1993), as well as 
the expectation that the ejection speed is proportional to 
the latitudinal variation of the escape speed on the sur- 
face of a rotating star, as noted by several investigators. 
By adapting these ideas to an explosion with continuum 
driving, however, our model does not suffer from the diffi- 
cult problem of WCD inhibition due to effects associated 
with non-radial forces in line-driven winds (Owocki et al. 
1998). Whether or not our model is apphcable depends 
on the nature of the episodic mass-loss event in any in- 
dividual case; the mass loss must be strong enough that 
it is not dominated by a line-driven wind during an out- 
burst (see Smith & Owocki 2006; Owocki et al. 2004). 

Previous models to explain the shapes of bipolar neb- 
ulae around rj Carinae and other LBVs differ substan- 
tially from ours. Most approaches have used hydrody- 
namic simulations of interacting winds, where a fast wind 
sweeps into an equatorial density enhancement to pro- 
duce the bipolar shape, or variations of that scenario 
(Frank et al. 1995, 1998; Dwarkadas & Balick 1998; 
Langer et al. 1999; Gonzalez et al. 2004a, 2004b). How- 
ever, none of these produced an equatorial disk with the 
same apparent age as the polar lobes. Soker (2004; and 
references therein) has discussed a complex model where 
accretion onto a companion star drives jets that shape 
the bipolar lobes, much as in similar models for plan- 
etary nebulae, but this model also fails to account for 
the equatorial disk. Other approaches involving non- 
spherical stellar winds are closer to our own, where the 
bipolar shape is an intrinsic feature of ejection from a ro- 
tating star (Owocki 2003, 2005; Owocki & Gayley 1997; 
Dwarkadas & Owocki 2002; Maeder & Desjacques 2001). 
In these models, the shape of the Homunculus is achieved 
by such a wind blowing with an enhanced mass-loss rate 
for a short time. Once again, however, none of these 
produced an equatorial disk.^ 

^ Maeder & Desjacques (2001) presented a second case where 
a steady stellar wind included a dense disk. However, while this 
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Fig. 4. — The initial kick velocity function (left panel), in units of ^GM / Rp, tuned to reproduce the shape of the Homunculus (center 
panel; solid and dotted lines have the same meaning as in Fig. [2J. A 3 dimensional rendering of the corresponding surface of revolution is 
shown in the right panel; the inclination (i = 41°) and orientation (P. A. = 310°) of the surface are based on the values published by Smith 
(2006). 



One model that did simultaneously produce bipolar 
lobes and a disk was the thermally-driven magnetohy- 
drodynamic rotating stellar wind model of Matt & Bal- 
ick (2004). While that model was encouraging, we felt it 
was also useful to pursue alternatives that did not rely 
on strong magnetic fields, since it is unclear whether the 
huge fields required to shape the 10-15 Mq (Smith et al. 
2003b) ejected during the Great Eruption are achievable. 

Two main assumptions required to produce the ob- 
served shape are explosive mass loss with continuum ra- 
diative driving and near-critical rotation. Explosive mass 
loss is justified by the observed fact that ry Car lost a 
huge amount of mass in a short time, with a required 
mass-loss rate that is too high to be accounted for with 
a line-driven driven wind, as noted earlier. Some recent 
stellar evolution models for very massive stars (Arnett 
et al. 2005; Young 2005) predict deep-seated hydrody- 
namic explosions that can potentially release the required 
amount of kinetic energy of almost 10^° ergs seen in the 
Homunculus (Smith et al. 2003b). Near critical rotation 
is expected to occur at late evolutionary stages as a BSG 
(e.g., Langer et al. 1999; Eriguchi et al. 1992). Soker 
(2004) has criticized single-star models for the ejection 
and shaping of the Homunculus based on the idea that 
the star spun down during the eruption and was not ro- 
tating sufhciently rapidly. However, two points of clari- 
fication should be mentioned in this regard. First, near 
the classical Eddington limit, a star can be near critical 
rotation for mass loss without necessarily being a true 
rapid rotator because of the effects of radiation pressure. 
Second, while the star may have shed angular momentum 
by ejecting the Homunculus and therefore was rotating 
more slowly after the outburst, it is the initial rotation 
and available angular momentum at the time of ejection 
that is critical, not the end state, especially in an explo- 
sion scenario such as that we are discussing here. 

The simultaneous formation of a disk and bipolar lobes 
from the same mechanism is a novel feature of our model, 
and there are interesting implications for the observed 

model may account for the enhanced density at the equator, it does 
not account for the shape (i.e. the speed) that can produce a disk 
at the same time as the polar lobes. It also depends on effects 
in line-driven winds that are not applicable to the Great Eruption 
of ?7 Carinae that produced the Homunculus; S. Owocki (private 
comm.) has noted the difficulty in forming disks via this type of 
opacity mechanism in radiatively-driven winds. 



structures. The equatorial ejecta of r/ Carinae were orig- 
inally described as a "skirt" rather than a "disk" because 
of their ragged, streaked appearance, which is different 
from the clumpy character of the polar lobes (e.g., Morse 
et al. 1998; Duschl et al. 1995; Zethson et al. 1999). 
This expanding excretion structure should not be con- 
fused with a Keplerian disk. In our proposed scenario 
for the formation of rj Car's skirt, material ejected from 
mid latitudes on the star migrates toward the equatorial 
plane, where it crashes into material from the opposite 
hemisphere. One can easily envision a violent collision 
that could lead to the ragged spray of ejecta seen today, 
rather than a smooth disk, and hydrodynamic simula- 
tions of this would be interesting. It seems likely that 
this violent splashing at the equatorial plane could create 
even faster disk material than the speeds seen in our sim- 
ulations where we simply canceled-out the vertical com- 
ponent of the velocity. To be sure, there are additional 
complexities in rj Car's equatorial ejecta that we do not 
even attempt to treat here. Most notable among them 
are the origin of the NN jet (Meaburn et al. 1993; Wal- 
born & Blanco 1988) and the presence of younger ejecta 
from the 1890 eruption that appear to co-exist with older 
ejecta in the skirt (Smith & Gehrz 1998; Davidson et al. 
2001; Smith et al. 2004b). 

Finally, we expect that our proposed scenario could 
have wider applications beyond the Homunculus around 
rj Carinae. A simple mechanism for how a rotating hot 
star can simultaneously produce equatorial and polar 
ejecta might be relevant for the famous triple-ring neb- 
ula around SN1987A (Burrows et al. 1995), providing a 
possible way to circumvent difficulties in explaining this 
object via the GISW scenario (see Smith 2007). Like- 
wise, it may apply to bipolar ejecta and rings around 
other blue supergiants in our galaxy such as Sher 25 
(Brandner et al. 1997) and the LBV candidate HD168625 
(Smith 2007), as well as LBVs in general. If extragalac- 
tic analogs of rj Carinae behave similarly in their out- 
bursts, then one might expect the so-called "supernova 
impostors" to be significantly polarized. The mechanism 
may also be applicable to the short lived emitting disks 
around B[e] stars (Zickgraf et al. 1986, 1996), or possibly 
other hot stars where episodic mass ejection is important. 

An encouraging property of our model is that in cases 
where near-Eddington radiative driving is less influential 
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than in rj Carinae (such as in lower-luminosity plane- 
tary nebulae), the mechanism proposed here generalizes 
to a situation that reproduces a simple bipolar pinched- 
waist morphology without an obvious fast disk. Thus, 
scenarios like those in the left column of Figure 2 may 
have wide application to observed morphologies of bipo- 
lar planetary nebulae, which generally lack such disks 
(Balick & Frank 2002). 

Some supernova explosions are seen to be intrinsically 
asymmetric; the bipolar supernova ejecta in SN1987A 
are seen directly (Wang et al. 2002), while others show 
polarization at early times (e.g., Leonard et al. 2000, 
2001; Leonard & Filippenko 2001). If near-critical rota- 
tion is important at some point within these explosions, 
it is conceivable that our model described here may of- 
fer an alternative to jet-driven hydrodynamic models for 
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explaining some aspects of asymmetric core-collapse su- 
pernovae. 

In any case, a viable mechanism for a single star to pro- 
duce strongly axisymmetric ejecta means that the pres- 
ence of asymmetry in the circumstellar environment is, 
by itself, not a valid justification to invoke close binary 
interac-tions in a supernova progenitor or any other hot 
massive star. 
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